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ABSTRACT 

We present the diffuse X-ray emission identified in Chandra observations of the young, massive 
Galactic star cluster Wester fund 1. After removing point-like X-ray sources down to a completeness 
Hmit of «2 X 10^^ erg s~-^, we identify (3± 1) x 10^'' erg s~^ (2-8 keV) of diffuse emission. The spatial 
distribution of the emission can be described as a slightly-elliptical Lorentzian core with a half-width 
half-maximum along the major axis of 25"±1", similar to the distribution of point sources in the 
cluster, plus a 5' halo of extended emission. The spectrum of the diffuse emission is dominated by 
a hard continuum component that can be described as a kT>3 keV thermal plasma that has a low 
iron abundance (^0.3 solar), or as non-thermal emission that could be stellar light that is inverse- 
Compton scattered by MeV electrons. Only 5% of the flux is produced by a kT^O.7 keV plasma. 
The low luminosity of the thermal emission and the lack of a 6.7 keV iron line suggests that <40,000 
unresolved stars with masses between 0.3 and 2Mq are present in the cluster. Moreover, the flux in 
the diffuse emission is a factor of two lower than would be expected from a supersonically-expanding 
cluster wind, and there is no evidence for thermal remnants produced by supernovae. Less than 10~^ 
of the mechanical luminosity of the cluster is dissipated as 2-8 keV X-rays, leaving a large amount of 
energy that either is radiated at other wavelengths, is dissipated beyond the bounds of our image, or 
escapes into the intergalactic medium. 

Subject headings: X-rays: stars, ISM — stars: winds — supernova remnants — star clusters: individual 
(Westerlund 1) 



1. INTRODUCTION 

Sensitive X-ray observations are an increasingly- 
important tool for studying young star clusters, partic- 
ularly now that the Chandra X-ray Observatory and the 
Newton X-ray Multi-Mirror Mission have made harder 
X-rays (2-10 keV) available for study. Young stars of all 
types are strong X-ray sources, with low-mass {M<3Mq) 
pre- main- sequence stars producing X-rays in their 
active magnetic cor onae fFrcibisch & Feigelson 2005; 
iFeieelson et al.ll2005j) . and massive OB stars {M>8Mq) 
producing X-ray s throug h shoc k s in their stellar winds 
ilChlebowski fc G armanvl I199H iBerghofer et alJ Il997t 
iSkinner et alTmoa) . Therefore, using observations of lo- 
cal star forming regions (e.g., Orion) as templates, mea- 
surements of the integrated X-ray luminosities of more 
distant clusters can be used to constrain their total stel- 
lar population, including the numbers of young stars that 
may be unobservable in the optical and infrared because 

of extinction or sour ce confusion (,Fcigclson et al 2005: 
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X-ray observations of clusters of massive stars also 
reveal diffuse X-ray emission that is produced as stel- 
lar winds encounter each other and the s urround- 
ing in terstellar mediu m (ISM: e.g.. iStevens fc Hartwe3 
12003"; 'Townslev et al.' "2003^ iLaw fc Yusef-Zadeh i200l 
Townslcv ct al. 2005b, 2006). Learning the fate of the 
energy carried by these winds, and eventually by super- 
novae, would provide insight into how galaxies evolve. 
If the energy is transferred to the ISM, it might at first 
trigger future generations of star formation, but a suffi- 
ciently large input of energy could clear away the ISM 
and halt star formation. Alternatively, if the energy es- 
capes a galaxy, stellar winds and supernovae would en- 
rich the intergalactic medium with metals. To determine 
the fate of that energy, it is necessary to obtain X-ray 
observations of clusters that have a range of ages and 
populations, and that are surrounded by ISM with a va- 
riety of densities fe.g.. iTownslev et al.> 2003. .2005b') . 

In this paper, we report on Chandra observations 
of the diffuse X-ray emission from the young Galac- 
tic star cluster Westerlund 1. The cluster con- 
tains 24 Wolf-Rayet (WR) stars, more than 80 blue 
super-giants, at least 3 red super-giants, a luminous 
blue variable, and an amazing 6 yellow hyper-giants, 
only 6 of which arc known in the ent ire rest of the 
Galaxy dW cstcrlund 1987: Clark & N egueruelal I2002L 
[2OOI iNemleruela fc Clark 2005: .Clark et al." "20051 As- 
suming a standard initial mass function lKrouDa,2002,) . 
Westerlund 1 could be as massive as 10^ Mq, mak- 
ing it several times larger than the well-known, young 
Galactic clusters the Arches, Quintuplet, and NGC 
3603. Westerlund 1 is also located only «5 kpc away 
(jClark fc Negueruelal 12002: Clark et a l.'2005). so it is 
one of the closest young, dense star clusters. Therefore, 
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Fig. 1. — Images of the 10' by 10' field around the center of Westerlund 1. The left panel displays the unbinned image, with the point-like 
X-ray sources described in Clark et al. (in prep.) marked with circles. The right panel displays an image of the diffuse X-ray emission in 
which the point sources have been excised and image has been adaptively binned to a signal-to-n oise of 10 using the Weighted Voronoi 
Tessellation algorithm i Dichl fc Statlcr 2006:1 . The location of the 10.6 s X-ray pulsar IMuno et alJl200(D is also indicated. 



Westerlund 1 is a crucial object for understanding the 
evolution of star clusters, and their impact on the ISM 
of their host galaxies. 

This is one of sev eral papers describ ing the Chan- 
dra observations. In iMuno et alJ l)2006() . we reported 
the detection of a s low X-ray pulsar in Westerlund 1. 
iSkinner et alJ 1)20061) examined the X-ray emission from 
the Wolf-Rayet (WR) stars in the cluster, as well as a 
subset of the OB supergiants that are brightest in X-rays. 
In Clark et al. (in prep.), we will report the spectroscopic 
identification of further optical counterparts to the X-ray 
sources, and discuss the origin of the X-ray emission from 
these stars. In this paper, we describe the spatial distri- 
bution f i)2.1|l and spectrum f H2.2|) of the diffuse X-rays. 
We compare the emission seen from Westerlund 1 to that 
of other massive star clusters in the Local Group (§3). 
We suggest that Westerlund 1 is one of only a few star 
clusters to produce mostly non-thermal X-rays. We dis- 
cuss the contraints that this places on the contributions 
of pre-main-sequence stars (§3.1), stellar winds (§3.2), 
and supernovae (§3.3) to the diffuse emission, and exam- 
ine what could be causing non-thermal emission (§3.4). 

2. OBSERVATIONS AND DATA ANALYSIS 

Westerlund 1 was observed with the Chandra X-ray 
Observatory A dvanced CCD Spectrom eter Spectroscopic 
array (ACIS-S; IWeisskopf etld]l200l on two occasions: 
on 2005 May 22 for 18 ks (sequence 5411), and on 2005 
June 20 for 42 ks (sequence 6283). We reduced the obser- 
vation using standard tools that are part of CIAO version 
3.3. We first created a composite event list for each ob- 
servation. We corrected the pulse heights of the events to 
mitigate for the position-dependent charge-transfer inef- 
ficiency using the standard CIAO process and excluded 



events that did not pass the standard ASCA grade fil- 
ters and Chandra X-ray center (CXC) good-time filters. 
We searched for intervals during which the background 
rate flared to > 3tT above the mean level, and removed 
one such interval lasting 3.6 ks from sequence 5411. The 
composite image of the full field is displayed in Clark 
et al. (in prep.); Figure displays the inner 10' of the 
cluster at 1" resolution. 

As described in Clark et al. (in prep.), we identi- 
fied point-like X-ray sources in each ob servation using 
a wa velet-based algorithm, wavdetect l|Freeman et alJ 
1200^ . In order to examine the diffuse X-ray emission, we 
then removed events that fell within circles circumscrib- 
ing approximately 92% of the PSF at the location of each 
point source, and created an image using the remaining 
photons. Within 5' of the cluster core, 7386 counts were 
associated with known point sources, and 38,350 counts 
in diffuse emission, so photons from point sources in the 
wings of the point spread function contribute only 0.5% 
to the diffuse flux. 

2.1. Spatial Distribution 

The signal-to-noise in a 1" pixel was low, so we 
adaptively binned the image using the Weighted 
Voronoi Tessellation algorithm implemented by 
Diehl fc St.atler (2 006), which is based on the algo- 
rithm of lCappellari fc Copiij l|2063i) . The resulting 
image is displayed in the right panel of Figure Q 

In order to quantify the extent of the diffuse emis- 
sion, we modeled its adaptively-binned, two-dimensional 
spatial distribution (Figure 1) as a Lorentzian function. 
Other functional forms used to model the light from 
star clusters also may be consistent with the data (e.g., 
lElson. Fall, fc Freem'anlll987l:IAnders. Gieles. fc de Griid 
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Fig. 2. — The radial profile of the diffuse X-ray emission (his- 
togram with errors) along with the best-fit model described by 
Equation |^ 



TABLE 1 
Lorentzian Model of the 
Distribution of Diffuse X-rays 



Parameter 



Value 



C 10-^ ph cm-2 s-i 

N 10-8 ph cm-2 s-i 

Qo (J2000) 

So (J2000) 

ro (arcsec) 

e 

8 (degrees) 



0.53±0.04 
3.4±0.2 
16 47 4.3±0.1 
-45 50 59±1 
25±1 
0.75±0.02 

13±3 
360/120 



Note. — We define 
rotations east of north. 



as positive for 



120061) . However, there is little tradition in modeling the 
spatial distribution of diffuse X-rays from star clusters 
with analytic functions, because that emission usuall y 
has a complex morphology fe.g.. iTownslev et alJl2003|) . 
so there is not an obvious choice for a functional form. 
We chose a Lorentzian function for its simplicity, and 
because it is similar to the King models often used to 
quan tify the di stribution of optical light from star clus- 
ters l|Kintdll962 '). The diffuse emission from Westerlund 
1 is not circularly symmetric, so we allowed for an ellip- 
tical distribution defined as 



/(x',y') = C + 



N 



1 + + eV')Ar 

[a — ao) cos 6o cos 9 + (6 — Sq) sin t 
y' = {5 — 6o) cos 9 ~ (a — ao) cos 5q sin I 



(1) 



X = {a — 



Here, ao and 5q are the center of the distribution and 
x' and y' are the offset in arcseconds from the center, 
where the axes defining them have been rotated east of 
north by 9 degrees. The remaining parameters are the 
background count rate C, the peak count rate TV, the 
ellipticity of the distribution e (a value of 1 implies a 
circle), and the characteristic radius of the distribution 
tq. The final parameters and the goodness-of-fit x^/v 
are listed in Table |21 

The fit is formally poor, with x^/i'^B, because there 
are structures in the diffuse emission that can not be de- 



scribed as part of an elliptical Lorentzian, including some 
bright knots of emission at the center of the cluster, and a 
ridge of emission extendin g to the so utheast toward the 
X-ray pulsar reported bv iMuno et al. (2006, it is also 
labeled in Fig. However, the model does provide a 
useful means to quantify the azimuthally-averaged dis- 
tribution of the emission, as can be seen in the radial dis- 
tribution, plotted in units of r' = (x'^ -|-e^y'^)^/^, in Fig- 
ure [3 The half- width half- maximum of the distribution 
is 25±1", which for a distance of 5 kpc corresponds to 0.5 
pc. The widths of the distributions of both optical stars 
and p oint-like X-r ay so urces are also «25" ijClark et alJ 
I2005t IMuno et al]l2006() . Moreover, the centroid of the 
diffuse emission lies within the 5" uncertainty in the cen- 
troid of the point sources, ijMuno et alJl200q) . although 
it is »20" from the cetroid of the optically-d etected stars 
l|Piatti. Bica. fc ClariallOOllClark et alfc oOS). The dis- 
crepancy could be caused either by differential extinction 
toward the cluster, or by substructure in the cluster (J. 
S. Clark et al., in prep.). 

In Figure [3 we have also indicated the amount of flux 
expected from the background of particles impacting the 
detector (1.5 x 10~^ ph cm"^ s~^), and the mean flux 
taken from observations of the Galactic plane at I — 28° 
and b = 0.2° (2.2 x 10"^ ph cm"^ s~\ o f which hah was 
particle background; see lEbisawa et all 12005). Even 5' 
from the cluster core, 60% of the flux was from Wester- 
lund 1, so there is a broad halo of diffuse X-rays around 
the cluster. In contrast, although a few O stars that are 
cluster members are located beyond ^3' from the clus- 
ter core ijClark et al„ 2005) . the surface density of X-ray 
point so urces beyond 3' is consistent that of the Galac- 
tic disk l|Muno et al.ll200^ . so there is no similar halo of 
point-like X-ray sources. 

2.2. Spectra 

Guided by Figure El we extracted spectra, response 
functions, and effective area curves for a circular 1' region 
around the cluster center, and annular regions 1-2', 2- 
3.5', and 3.5-5' from the cluster center. The background- 
subtracted spectra are displayed in Figure |3| in units of 
detector counts per square arcminute. 

The spectra contained contributions from at least three 
sources of X-rays: plasma and unresolved stars in the 
cluster, diffuse emission from the Galactic plane, and 
events produced by particles incident on the detectors. 
The spectrum of the background from particles has been 
well-characterized using observations in which the ACIS 
detectors were stowed out of the focal plane. Therefore, 
we subtracted the spectrum of the particle background 
from our source spectra. However, we were not able to 
make a local estimate of the Galactic emission, because 
the cluster emission extended over the entire image. We 
have not attempted to subtract the Galactic plane emis- 
sion from the spectra of the diffuse emission from Wester- 
lund 1, but have estimated the contribution of the Galac- 
tic emission by modeling a spectru m from observations 
of a field at I = 28° and b = 0.2° l|Ebisawa et al.ll20011 

We modeled the emission using version 12.2.0 of 
XSPEC (Arnaud et al. 1996). We chose to model only 
the 1.5-8 keV energy range, for two reasons. First, the 
mean absorption column measured from the X-ray spec- 
tra of the point sources in Westerlund 1 is equivalent to 
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Fig. 3. — Background subtracted spectra of the diffuse X-ray emission from Westerlund 1. In the top fialf of each panel, the spectra are 
plotted in detector units, so the intrinsic spectrum is still convolved with the response of the detector and telescope. The model spectrum 
is plotted with the solid histogram. The dotted line denotes the thermal component of the model, and the dashed line the non-thermal 
power-law component. The spectrum becomes noticeably harder farther from the cluster, because the thermal component becomes weaker. 
The grey data points represent Galactic plane emission at Z = 28° and b = 0.2°, which contributes only 30% of the total flux toward 
Westerlund 1. In the bottom half of each panel, we plot the difference between the data and model normalized to the uncertainty (cr) in 
the data. 



2 X 10^^ cm^^ of hydrogen (J. S. Clark et al., in prep.). 
Therefore, most of the observed X-ray flux from the clus- 
ter should emerge at energies above ~2 keV, and the 
lower-energy X-rays are probably from foreground emis- 
sion. Second, given that the diffuse emission from West- 
erlund 1 is probably from several different sources, we 
do not have enough physical guidance to extrapolate our 
models below 2 keV. When we try to apply simple mod- 
els, the inferred de-reddened 0.5-2.0 keV flux can span an 
order-of-magnitude depending upon the model assump- 
tions. The lower bound was chosen to include the promi- 
nent line at 1.8 keV in the spectra from the inner 2' of 
the cluster. For consistency with the other works quoted 
in §3, we report the observed and de-reddened 2-8 keV 
fluxes. 

We first attempted to model the spectrum as a sin- 
gle temperature thermal plasma, either in or out of coUi- 
sional ionization equilibrium, absorbed by the interstellar 
medium. We found that those simple models provided 
a poor description of the data from the inner 2' of the 



cluster (x^/j^ > 1.5). 

Therefore, we modeled the spectra as an absorbed, 
two-temperature thermal plasma. Most of the 4-8 keV 
continuum flux could be modeled as a hot plasma with 
kT2^ 3 keV. In the inner 3.5' of the cluster, the pres- 
ence of emission lines near 1.8 keV from He-like Si and 
2.3 keV from He-like S indicated that some of the flux 
is produced by a cool kTi<l keV plasma. The metal 
abundances in the cooler component were poorly con- 
strained because it contributes very little to the con- 
tinuum emission, so we fixed the metal abundances to 
the mean best-fit value of Z/Zq—2. Moreover, the spec- 
trum taken from the 3.5-5.0' annulus lacked obvious line 
emission, so the parameters of any cool plasma emis- 
sion were unconstrained. Therefore, we omitted the cool 
component from the model of that spectrum. In Ta- 
ble El we list the parameters of the best-fit, two tem- 
perature, coUisional ionization equilibrium models; us- 
ing non-equilibrium models yields similar results for the 
derived temperatures and abundances. Using these as- 
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TABLE 2 

Two-Temperature Plasma Model for the Diffuse X-ray Emission 



Parameter <1' 1-2' 2.0-3.5' 3.5-5.0' 



counts 4259 6679 12698 14714 

background 698 2244 5613 6636 

area (arcmin^) 2.7 8.9 25.6 34.9 

Nn (1022 cm-2) 2.2t°-^ 2.0^^2 2.2t°-^ 2.4t°-} 

fcTi (keV) 0.71 0.7tli l-OtH 

Zi/Zq 2.0 ' 2.0 ' 2.0 ' 

/ rie iriH idV (cm-6 pc^) 9^1 ^tl ^tl 

fcT2(keV) 3.2tO-j 5.7tJ;2 11^^ 6.3«;«g 

Z2/Z0 <0.4 <0.3 0.6t°j <0.2 

/ n.e,2n/j-,2dV (cm-6 pc^) 16j;2 17^2 241:2 4313 

X^/i^ 73.4/65 80.4/72 104.9/101 129.1/116 

Fx (10-^3 erg cm-2 s-i) 7.1±0.9 11.3±1.0 23.8±2.0 27.0±2.0 

M^x.i (10-1^ erg cm-2 s-i) 1.5 0.9 2.0 

ttFx'2 (10~^^ erg cm-2 s"^) 8.3 1A4 27^7 34.9 

Note. — Uncertainties are Icr, found by varying each parameter until Ax^ = 
1.0. ui^x.l and 'ui^x,2 are the de-absorbed 2-8 keV flux from the thermal and 
non-thermal components of the spectral model, respectively. If we extrapolate 
our models to the 0.5-2.0 keV band, the rapidly-increasing contribution from 
the fcT=0.7 keV thermal plasma causes the inferred X-ray luminosity to be a 
factor of 2-3 larger. Note also that at 5 kpc, 1 arcmin = 1.45 pc. 



TABLE 3 

Thermal Plus Non-Thermal Model for the Diffuse X-ray Emission 



Parameter <l' 1-2' 2.0-3.5' 3.5-5.0' 



TVh (1022 cm-2) 2.71°;^ 2.3j;»;^ 2.4j:E;;2 2.8+°-^ 



0.2 ^-^-0.2 ^■°-0.2 

-I-O-IO n «a+0-09 1 1+0.2 



fcTi (keV) 0.68+^;i^ O.mZrvi l-l-o.i 

Zi/Zq 2.0 2.0 2.0 

/nennrfV (cm-6 pc3) lll^ 7.ll? * 6^2 

P 97+0.2 91+0:1 17+0.1 9 n+O-2 

^ ^-'-0.2 ^-^-0.2 ^-'-0.1 ^■'-'-0.1 

iVr (10-3 ph cm-2 s-i ^eV-i) 1.01°^ 0.7lE!:2 0.sto\ 1.7t°;* 

X^/v 70.4/66 80.4/73 107.2/102 132.2/116 

Fx (10-13 cm-2 g-i) 7.2±0.7 11.7±1.0 23.1±2.0 28.0±2.0 

M-Fx,i (10-13 erg cm-2 s-i) ^ q 2.2 

mJx!2 (10-" erg cm-2 g-i) g^g 1A2 27^2 37.3 

Note. — See Table 2. 



sumptions, the models were generally good descriptions 
of the data, with x^/vKil. 

The most notable trend in the cool components is that 
their contributions to the spectra decline from 15% in the 
central 1', to 7% between 1' and 3.5', finally becoming 
undetectable beyond 3.5' from the cluster center. Other- 
wise, the inferred interstellar absorption remains roughly 
constant near 2.6 x 10^^ cm"^, the temperature of the 
thermal component is constant near fcT«0.7-l keV. We 
find that the temperature of the hot plasma increases 
from kT2—3 keV at the cluster center, to a maximum 
11 keV in the 2.0-3.5' annulus, and then decreases to 
6 keV in the outer annulus. The relative lack of flux 
near the He-like Fe line at 6.7 keV in most of the spec- 
tra implies that the iron abundances are less than half 
of the solar values. Interestingly, similar sub-solar iron 
abundances are inferred from the lack of 0.8-1.0 keV Fe 
L lines from several known Galactic Wolf-Rayet and O 



stars (e.g.. ISkinner et alJl200ll 12001 ISchulz et alJl2003l: 
ISkinner et a l."'2005'). 

Alternatively, the lack of lack of line emission near 6.7 
keV could be explained if much of the continuum X-ray 
emission is non-thermal. Therefore, we have also mod- 
eled the emission as the sum of emission from a kT<l 
keV thermal plasma and a power law (Table 13) . The 
metal abundances in the cooler component were once 
again fixed to Z/Zq—2, and we omitted the cool com- 
ponent from the model of the 3.5-5' annulus. This pro- 
vides an equally good description of the data as the two- 
temperature plasma model, and the same trends are ev- 
ident: the contribution of the cool component declines 
monotonically with offset from the cluster center, and 
the overall spectrum becomes harder. 

The contributions of each model component under the 
second set of models are indicated in Figure^ using dot- 
ted lines for the thermal plasma, and dashed lines for the 
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power law component. We also display the spectrum of 
the Galactic ridge emission at I = 28° and b — 0.2° (grey 
data points). The line emission from the Galactic flux is 
a bit stronger than that from Wester lund 1, but other- 
wise the spectra are fairly similar. Therefore, we cannot 
completely rule out the hypothesis that the emission be- 
yond «2' from the cluster core is Galactic. However, our 
assumption that the diffuse emission is from Westerlund 
1 is conservative. As described in §3, we find that the 
luminosity of diffuse X-rays from Westerlund 1 is much 
lower than expected, and assuming that the diffuse halo 
is Galactic would exacerbate the discrepancy. 

The total, de-reddened 2-8 keV flux from within 5' of 
Westerlund 1 is 9.3 x 10~^^ erg cm^^ s^^. By varying 
the assumptions in our model, we find that the system- 
atic uncertainty in the 2-8 keV flux is «20%. Based on 
the Chandra o bservations take n at I = 28° and b = 0.2° 
(see also .Ebisawa et all l2005j) , we expect the Galactic 
emission to be 3 x 10~^^ erg cm~^ s~^ arcmin"^ (2-8 
keV; see also Hands et al. 2004). Therefore, within 5' of 
the core of Westerlund 1 the Galactic plane contributes 
«20% to the inferred flux. Subtracting this foreground 
and background emission, and using a distance to West- 
erlund 1 of 5 kpc (Clark ct al. 2005), we find that the 
luminosity of the diffuse X-ray emission from the cluster 
is (3 ± f) X 10^4 erg s"! (2-8 keV). Only «5% of this 
luminosity is from the <1 kcV thermal component. 

3. DISCUSSION 

The origin of the diffuse X-ray emission from clus- 
ters of massive young stars is currently under debate. 
Several authors (e.g., Ca nto, Raga, & Rodriguez 2000; 
iStevens fc HartwelM 12003(1 have modeled the diffuse X- 
rays from the most massive clusters in the Local Group 
as a cluster wind. Under this model, the winds of individ- 
ual stars collide, thermalize, and form a pressure-driven 
bulk flow that expands supersonica ll y into the interstellar 
mediu m ((Chevalier fc Clegel 11985(1 . IStevens ^ HartweU 
((2003() tabulated results from studi es of R13 6, NGC 
3603 f Moffat et al."200 a), NGC 346 ((Naze et al.ll200l . 
the Rosette ( Townsl ey et alJ l2003() . and the Arches 
(lYusef-ZadelTet al.l (200^ and showed that the lumi- 
nosities of their diffuse X-ray emission ([1 — 6] x 10'^^ 
erg s~^) were considerably larger than would be ex- 
pected from the standard cluster wind model. The large 
X-ray luminosities can be explained several ways: the 
densities of the cluster winds could be higher than ex- 
pected because the stellar winds entrained cooler mate- 
rial or because radiative los ses decreased the tempera - 
ture of the shocked plasma ((Stevens Sz HartweU l2003() : 
the wind energy could be dissipated through heat con- 
duction where it encounter s nearby molecular clouds 
((Dorland &: Montmerlelir987(l : or the wind could be con- 
fined by the surrounding ISM ((Chu et alJll995|) . Alter- 
natively, the large X-ray luminosities might partly result 
from the fact that unresolved pre-main-sequence stars 
should contribute significantly to the luminosity of the 
(apparently) diffuse emission, espec ially for more distant 
clusters (^e.g..lTownslev et alJl200^. 

Westerlund 1 is a t lea st as ma ssive as NGC 3603, R136, 
and the Arches (Clark et al. 2 0051) . so from an obser- 
vational standpoint the luminosity of its diffuse X-rays 
([3±1] X 10^" erg s^^; 2-8 keV) is understandable. How- 
ever, the spectrum and spatial distribution of the dif- 



fuse X-ray emission from Westerlund 1 presents more of 
a puzzle. First, the spectrum lacks the line emission 
from He-like Fe that would be expected from a ther- 
mal plasma given the hard continuum flux. In contrast, 
the spectra of pre-main sequence stars exhibit prominent 
lines from He-like and H-like Si, S, Ar, and Fe that im- 
ply n ietal abundances up to ten times the solar value 
('e.g.. iFeigelson et"al]l2005() . However, X-ray spectra of 
O and WR stars often exhibit we ak Fe emission that 
imply abimdances <0.3 sol ar (e.g., [S kinner e t alJl200lL 
2002; Schulz et al. 2003; Skinner et al . _20Q5), so it is 
possible that the diffuse flux is dominated by plasma 
from the O and WR star winds, with little contribu- 
tion from pre-main-sequence stars. It is also possi- 
ble that the diffuse emission is non-thermal, by anal- 
ogy with similar interpretations for the hard flux from 
a h andful of young s tellar associations, in cluding RCW 
38 (IWolk et alJl2(ffll. DEM L192 rN 51D: ICooner et all 
I2004() . 30 Dor C llBamba et alJl2004fl. and possibly the 
Arches cluster ((Law fc Yiisef-Zadehll2003) . In most of 
the above cases, the non-thermal emission has been in- 
terpreted as synchrotron emission from supernova rem- 
nants. We will examine this hypothesis for Westerlund 1 
in S(23| 

The second surprise is that the the diffuse X-ray emis- 
sion from Westerlund 1 seems to extend far beyond the 
core of the cluster. Within the inner 2', the surface 
brightness of the diffuse emission falls off with a half- 
width half-maximum of 0.5 pc (Fig. ^ and O, which is 
identical to the distribution of point-like X-ray sources 
(Clark et al., in prep). This core of diffuse emission could 
be produced either from the cluster wind, which radiates 
X-rays mostl y in the region wh ere the colliding winds are 
thermalized ((Stevens fc Hartw eU 2003)j O'^ from an unre- 
solved population of pre-main sequence stars. However, 
between 2' and 5' (3-7 pc) from the cluster core the dif- 
fuse X-ray flux attains a constant level x 10~^^ erg 
cm~^ s~^ arcmin"^ (Tab. [2 Fig. which is 2-3 times 
larger than is exp ected from the Galactic plane (e.g., 
iHands et all 120041: lEbisawa et ani200 5'l. An expanding 
thermal plasma would exhibit a rapidly-declining tem- 
perature profile, yet the spectrum of this halo of diffuse 
emission is quite hard and lacks the line emission ex- 
pected from a cooling plasma. This makes it tempting 
to interpret the diffuse halo as non-thermal particles that 
are accelerated in a large-scale outflow. 

Therefore, although the luminosity of the diffuse X-ray 
emission from Westerlund 1 is not surprising, the lack of 
line emission in the spectrum and broad spatial distri- 
bution of the diffuse X-rays is. To address this, in the 
following sections we quantify the probable contributions 
of pre-main-sequence stars, stellar winds, and supernovae 
to the X-ray emission from Westerlund 1. 

3.1. Unresolved Low-Mass Stars 

The non-thermal spectrum of the diffuse X-ray emis- 
sion from Westerlund 1 puts interesting constraints on 
the population of low-mass stars in the cluster. The aver- 
age spectrum of the lightly-absorbed pre-main-sequence 
stars in Orion can be qualitatively described as a two- 
temperature plasma with fcTi=0.5 keV and fcT2=3.3 keV, 
and with metal abund ances of up to 10 times solar for 
S, Ar, Ca, and Fe (,Feigelson et alJl2005D . In contrast, 
in our models for the diffuse emission from Westerlund 
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1, a kT<l keV plasma contributes only 5% of the 2-8 
keV diffuse X-ray flux, and the remaining hard flux does 
not exhibit the expected He-like Fe line at 6.7 keV, plac- 
ing an upper limit on the Fe abundance of <0.3 solar 
(Tab.EJ. To obtain a conservative estimate of the num- 
ber of low mass stars in the cluster, we assume that they 
have solar Fe abundances (i.e., much lower abundances 
than the stars in Orion). We find that implies that they 
they produce <30% of the diffuse flux from Westerlund 
1, or «9 X 10^3 erg s-^{2-8 keV). 

We use the results of the Chandra Orion Ultradeep 
Project (COUP) to convert this luminosity into a number 
of low-mass stars, taking into account the difference in 
ages between the two clusters. For Orion, the 1398 stars 
later than B4 in the COUP observations have an inte- 
grat ed, de-absorbed 2-8 k eV luminosity of 1.2 x lO'^'^ erg 
s-i llFeigelson et al.ll2005j) . Most of this emission is pro- 
duced by stars with 0.3<M<3 Mq.^ However, the star s 
in Westerlund 1, with ages of «4 Myr (|Clark et alJ2005fl . 
are significantly older than the 1-Myr-old population in 
Orion. To take this into account, we first note that when 
2-3M0 stars reach an age of ~4 Myr, they become fully 
radiative and their X-ray luminos ities drop by an order 
of magnitude ijFlaccomio et al.|l20 63) . Even though they 
are only 5% of low-mass stars by number, in Orion these 
2-3Mq stars produce «30% of the flux from 0.3<A/<3 
Mq stars (see Fig 4. in Feigelson et al. 2005), or ^8x10^^ 
erg s~^. Second. .Preibisch fc Feigelson 1,2005^1 find that 
the X-ray luminosities of young stars with 0.5<M<2Mq 
fall off with time t as Lx oc r""'^^, so at 4 Myr the stars 
in Orion should be ~3 times fainter. Therefore, if Orion 
were 4 Myr old, we would expect its «1400 stars with 
0.3<M<2 Mq to have a luminosity of 3 x 10^^ erg s"^ 
(2-8 keV). Our upper limit to the integrated X-ray lu- 
minosity of low-mass stars in Westerlund 1 is <9 x lO'^^ 
erg s~^, so we infer that Westerlund 1 contains <40,000 
stars with masses between 0.3<M<2Mq. 

This number of low-mass stars is smaller than one 
would expect if one were to extrapolate from the number 
of massive, post-main-sequence stars in the cluster using 
a standard initial mass function (Kroupa 2002). There 
are «150 stars brighter than V=21 within 5' of the cen- 
ter of Westerlund 1, the faintest of which have recently 
been identified as 07 main sequence stars that would 
have initial masses >30 Mq (J. S. Clark et al. in prep). 
The maximum initial mass of the stars remaining in the 
cluster is uncertain because there is no precise means to 
determine the initial m asses of the supergiants and WR 
stars, but IClark et"al] |2005) argue that it is probably 
in the range of 40-50 Mq . If we assume that the initial 
mass function can be described as a broken power law of 
the form dN oc M~"dM^, where a=2.3 f or M.,>0.5Mq 
and a=1.4 for 0.3<Mi<0.5 Mq llKroupat.2002,) , then if 
there are 150 stars with 30<Mi<50 Mq, we would ex- 
pect 100,000 stars with 0.3<M<2 Mq. This inferred 
lack of low-mass stars can be explained several ways. The 
slope of the initial mass function could be flat (a<2.1), 
as has been inferred for NGC 3603 and the Arches clus- 

® The COUP sample of X-ray sources is complete to !^0.2Mq, 
and excluding the singularly bright 06 star 9^ Ori C (with a lu- 
minosity of 3 X 10^^ erg s~^ [2—8 keV], it could be detected as a 
point source in our observations of Westerlund 1) the >3Mq stars 
in Orion produce only fa6% of the integrated X-ray luminosity. 



ter based on infrared sta r counts fe.g.. lEisenhauer et alJ 
119981 IStolte et al.ll2005j) . The mass function could be 
truncated at low masses {M<O.QMq), by analogy with 
the fact that M <7M(7) stars appea r to be depleted in the 
Arches cluster (fS tolte et al]l2005j) . Finally, the initial 
masses of the post-main sequence stars in Westerlund 
1 could span a much wider range o f masses (20-60Mo) 
than assumed in IClark 

If we assume that the mass function is truncated at 
low masses, the total mass of the cluster wo uld not differ 
signifi cantly from the estimate of 10^ Mq in IClark et alJ 
(j2005j) based on the un-modified Kroupa form. However, 
if the mass function is flat, or the optically-detected stars 
had a wider range of initial masses, the total mass of the 
cluster would be only 40,000-70,000 Mq. Obviously, an 
accurate measurement of the mass function, and conse- 
quently the total mass, requires direct infrared observa- 
tions of the low-mass stars in Westerlund 1. However, 
these X-ray observations provide a useful starting point. 

3.2. Stellar Winds 

It is not clear whether stellar winds or supernovae are 
the dominant source of the diffuse X-ray emission from 
Westerlund 1, because it is at an age when both should 
contribute equally to its mechanical output (e.g., Lei- 
therer, Robert, & Drissen 1992). Individual stellar winds 
carry a power of 

(2) 

where M is the mass loss rate, and Ww is the wind ve- 
locity. The WR stars dominate the mechanical output 
from stellar winds, with typical A/ « 6 x 10~^ Mq yr~^ 
and fw » 1700 km s~^, so that « 5 x 10^'^ erg s~^ 
ijLeitherer et al. "'l992l. With at 24 WR stars in the clus- 
ter, the total mechanical energy output from winds is 
> 1 X 10^9 erg s"^ 

We can estimate the X-ray luminosity of the resulting 
cluster wind using t he analytic solution s to the density 
and temperature that lCanto e t al. ( 2000| see also Stevens 
& Hartwell 2003) derived for a win d expanding super- 
sonic ally into the interstellar medium ijChevalier fc Clegel 
119851) . Within the radius of the cluster where the stars 
input their energy, the density is given by 

no = O.IA. ( f - , ] ( ' (^) " cm- 

' \^ lO-^Mo/yr J V km/s J \pc J 

(3) 

(note that we used the supersonic solution with an adi- 
abatic index 7 = 5/3 for the original equation) and the 
temperature by 

kTo = 1.3 f ^^^^^ keV, (4) 

° VlOOOkm/sy ' ^ ^ 

where N is the number of stars, M^ is the average 
mass loss rate of the stars, and Rc is the radius within 
which the stars are contained and the winds are ther- 
malizcd. If we use the values above for WR stars, as- 
sume iV=24, and take the radius of the cluster to be 
i?c=4 pc («3'at 5 kpc), then wc find no=0.6 cm^^ g^^if]^ 
kTo=4: keV. Using any of the standard plasma models 
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in XSPEC (e.g.. iMewe. Lemen fc van den OordI 119861: 
iLiedahl. Osterheld. fc GoldsteinI I1995D . and converting 
no and Rc to an emission measure (i.e., Kem = f^'^^c'^o) 
we find a predicted Lx = 3 x 10^"* erg s~^ (2-8 keV). 
Therefore, a cluster wind could in principle account for 
all of the diffuse X-rays from Westerlund 1 . 

However, the cluster wind model is not able to ac- 
count for the spatial distribution of the diffuse X-rays 
from Westerlund 1. A cluster wind would produce al- 
most all o f the diffuse X-ray emiss ion within the core 
radius Rc l|Stevens fc Hartwelill2003j) . whereas at 70% of 
the diffuse emission from Westerlund 1 is part of a halo 
that extends out to at least 5' (Fig. If consider only 
the core of the diffuse emission as originating from stellar 
winds, then the cluster is underluminous by a factor of 
two. 

This result is particularly surprising given that 
the s tandard cluster wind model applied by 
IStevens fc Hartwell (2003) to NGC 3603, R136, 
and NGC 346 predicted si gnificantly less flux t han is 
observed. In those cases. Tstevens fc Hartwelll l|2003D 
favored the hypothesis that cold material was being 
entrained in the wind. To reconcile our results for 
Westerlund 1 to the cluster wind model, we would 
have to assume that enough cold material is entrained 
that the cluster wind no longer emits in the 2-8 keV 
bandpass. This requires that the plasma be cooled 
by a f actor of wlO, to <0.4 keV. Based on Figure 1 
in ISte vens fc Hartwell (200;f), we estimate that this 
would require that the mass of cold material input into 
the wind is twice that of the hot material, or roughly 
3 X 10"^ Mq yT~^. The mass loss rates from the «10 
luminous blue variables, red supcrgiants, and yellow 
hypergiants, which should be <10^^ Mp, yr~ ^ each 
(e.g.. i.Tura fc K leinmann 1990; ILeitherer et alJ Il992t 
iSmith. Vink. fc de Koten .20041 . probably could not 
account for this large amount of cool mass. Therefore, 
either there is a currently-unseen source of mass in 
Westerlund 1, or the stellar winds are not thermalized 
within the cluster and escape without radiating much. 

3.3. Supernovae 

The presence of an isolated X-ray pulsar in Westerlund 
1 ijMuno et al.ll200^ confirms that supernovae have oc- 
curred there. If we extrapolate an initial mass function 
with slope a=1.8-2.7 for M>30Mq to higher masses, 
we expect that the cluster originally contained Ri80-150 
stars with initial masses >5OM0 that have already un- 
dergone supernova. For the most massive stars, this 
would have started when the cluster was about 3 Myr 
old, so the average supernova rate over the last 1 Myr 
should be on order one every 7,000-13,000 yr. If each 
supernova had a kinetic energy of 10^^ erg, then the av- 
erage power released by supernovae is ~ (2 — 5) x 10'^^ 
erg s~^. 

No obvious supernova remnant is present in our Chan- 
dra image of Westerlund 1, but this is not surprising. 
We have examined images from the Spitzer GLIMPSE 
program (R. Indebetouw, private communication), and 
there is no evidence that dense gas or dust still surrounds 
Westerlund 1. Therefore, Westerlund 1 appears to have 
cleared away the ISM for parsecs around. When a super- 
nova occurs in such an evacuated cavity, a typical radio 



and X-ray remnant is not expected until the remnant 
encounters t he boundaries of the bu bble blown by the 
cluster (e-g.- lCiotti fc D'EcroldlT989|) . 

Whether the hard, possibly non-thermal emission from 
Westerlund 1 is produced by supernovae is unclear. 
Most supernova remnants produce thermal X-ray emis- 
sion with strong lines, but a few are also non-thermal 
X-ray sources. For example, RCW 86 and SN 1006 
exhibit non-thermal filaments near the outer bound- 
ary of t he shock, and thermal emi ssion in the inte- 
rior ( e.g. jDver. Revnolds. fc Borkows ki 2004: R ho et all 
\mm . AX J1843.8-0352 and G346.3-0.5 exhibit non- 
thermal emissi on almo st exclusivel v throughout the 
remnant (e.g.. lUeno et al. 2003: La zendic et al.l l2005t 
iHiraga et all 12005(1 . Unfortunately, there is not a sat- 
isfactory explanation as to why a small fraction of su- 
pernova remnants produce non-thermal emission, so the 
issue remains unresolved for Westerlund 1. 

3.4. Non-Thermal Particles 

In principle, non-thermal parti cles can be produced ei- 
ther in supernova remnants (e.g., Lvutikov fc Pohl 2004) 
or in colliding stellar winds (e.g.. .Eichler fc Usov..l993j) . 
Once they are produced, inverse-Compton scattering 
should dominate synchrotron losses by a large factor in 
Westerlund 1 (see, e.g., Rvbicki fc Lightman 1979). The 
ratio of the energy-loss rates is given by the ratio of 
the background radiation to the magnetic energy den- 
sity. The energy density of the stellar light from the OB 
and WR stars in Westerlund 1 is approximately 

where Lgtars ~ IO^Lq the luminosity of the cluster. For 
synchrotron losses to be important, magnetic fields would 
have to have an energy of B"^ /Sir > ?7phot, which corre- 
sponds to B >0.4 mGauss. This is much stronger than 
the microGauss fields generally assumed for the interstel- 
lar medium (e.g.. Beck 2001), so inverse-Compton scat- 
tering is probably the dominant loss mechanism for non- 
thermal particles. 

If the non-thermal X-ray emission is produced by 
inverse-Compton scattering, the energy requirements are 
modest. Non-thermal particles would only need to be re- 
plenished at a rate sufficient to balance the X-ray lumi- 
nosity, 3 X 10'^'' erg s~^. Furthermore, inverse-Compton 
scattering photons from optical and UV energies {Ein = 
2-20 eV) into the X-ray band (Kut ~ 3 keV) only 
requires electrons with 7^ £^out/-E'in, or energies of 6- 
20 MeV. These particle energies are rather small. For 
comparison, if the magnetic field in Westerlund 1 has a 
strength of only 10 /iCauss, producing non-thermal syn- 
chrotron radio emission like that seen from 30 Dor C 
(Baniba ct al. 2004) requires electrons with energies of a 
few GeV. Therefore, detecting diffuse, non-thermal ra- 
dio emission from a star cluster like Westerlund 1 (e.g., 
lYusef-Zadeh et aLl l2003) would provide a much more in- 
teresting constraint on the maximum energies of parti- 
cles than detecting non-thermal X-rays. The int erfero- 
metric radio observations in the literature (Clark ~t alJ 
Il998|) would have resolved out arcminute-scalc diffuse ra- 
dio emission, so single-dish observations are necessary to 
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determine wliether higher-energy particles are also pro- 
duced by the cluster. 

4. SUMMARY 

We have identified diffuse X-ray emission within 5' of 
the core of Westerlund 1 with a modest luminosity of 
(3 ± 1) X 10^'* erg s"^ (2-8 keV). This low luminosity is 
puzzling, because unresolved pre-main-sequence stars, a 
thermalized cluster wind, or a series of supernova rem- 
nants would each be expected to produce at least this 
much X-ray emission. Therefore, one or all of these mech- 
anisms is not producing nearly as much X-ray flux as 
would be expected based on comparison with other star 
clusters and with theoretical calculations. 

The lack of a 6.7 keV He- like Fe line accompanying 
the hard 4-8 keV continuum implies that no more than 
30% of the diffuse emission is produced by young stellar 
objects. Therefore, we infer that there are <40,000 stars 
with masses between 0.3 and 2 Mq, which is significantly 
fewer than the 10^ stars one would expect from extrap- 
olating the number of massive, optically-identified stars 
to lower masses u sing a standard initial mass function 
ijClark et alJl2005l) . Moreover, this limit is conservative, 
because in computing it we have assumed that the line 
emission from the low-mass stars would be produced by a 
solar abundance of iron. If we had assumed that iron had 
an abundance several times t he solar value, as it d oes in 
the spectra of stars in Orion ijFeieelson et al.ll2'005() , then 
only a few percent of the diffuse 2-8 keV flux could be 
produced by pre-main sequence stars. 

In contrast, the lack of iron emission in the spectrum 
is consistent with a similar under-abundance of iron that 
is observe d in X-ray spectra of individ ual O and WR 
stars re.g.. lSkiTiTier et aljl200Tl [2002; Sc hulz et a,].ll20M 
iSkinner et all I2005H . However, if the O and WR star 



winds collide and thermalize as expected, they would 
form a pressure- driven cluster wind that would expand 
and c ool rapidly ijCanto et al.ll200(i iStevens fc HartwelH 
120031) . Such a wind would not radiate in the X-ray band 
outside of the cluster core, and therefore cannot explain 
the broad halo of emission between Ri3' and 5'. 

Instead, the halo of X-rays may represent non-thermal 
particles accelerated by the colliding stellar winds or by 
supernova remnants. However, the energy lost in X-rays 
represents less than 10~^ of the kinetic energy released 
by stellar winds and supernova remnants. The rest of 
the energy either (1) emerges below 2 keV where our 
observations are insensitive, (2) dissipates beyond the 
bounds of our image («7 from the cluster core) when the 
cluster wind or supernova remnants impact the ISM, or 
(3) escapes the Galactic plane, enriching the intergalactic 
medium with metals. We plan to address the second 
option in the next year, by observing a larger area around 
the cluster («15') with XMM-Newton. 
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